Chaos driven fusion enhancement factor at astrophysical energies 
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We perform molecular dynamics simulations of screening by bound target electrons in low energy 
nuclear reactions. Quantum effects corresponding to the Pauli and Heisenberg principle are enforced 
by constraints. We show that the enhancement of the average cross section and of its variance is 
due to the perturbations induced by the electrons. This gives a correlation between the maximum 
amplitudes of the inter-nuclear oscillational motion and the enhancement factor. It suggests that 
the chaotic behavior of the electronic motion affects the magnitude of the enhancement factor. 
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The knowledge of the bare nuclear reaction rates at 
low energies is essential not only for the understanding 
of various astrophysical nuclear problems, but also for 
assessing the effects of the host materials in low energy 
nuclear fusion reactions in matter. This is currently a 
subject of great interest in nuclear physics, since Muen- 
ster group has reported that the experimental cross sec- 
tions of the 3 He(d,p) 4 He and of 2 H( 3 He,p) 4 He reactions 
with gas targets show an increasing enhancement with 
decreasing bombarding energy with respect to the values 
obtained by extrapolating from the data at high ener- 
gies 0. Many studies attempted to attribute the en- 
hancement of the reaction rate to the screening effects 
by bound target electrons. In this context one often es- 
timates the screening potential as a constant decrease of 
the barrier height in the tunneling region through a fit 
to the data. A puzzle has been that the screening po- 
tential obtained by this procedure exceeds the value of 
the so called adiabatic limit, which is given by the differ- 
ence of the binding energies of the united atoms and of 
the target atom and it is theoretically thought to provide 
the maximum screening potential |2| . Over these several 
years, the redetermination of the bare cross sections has 
been proposed theoretically Q and experimentally Q, 
using the Trojan Horse Method The comparison 

between newly obtained bare cross sections, i.e., astro- 
physical S-factors, and the cross sections by the direct 
measurements gives a variety of values for the screening 
potential. There are already some theoretical studies per- 
formed using the time-dependent Hartree-Fock(TDHF) 
scheme HQ. 

In this letter we examine the subject within the con- 
strained molecular dynamics (CoMD) model 0, even in 
the very low incident energy region not reached experi- 
mentally yet. At such very low energies fluctuations are 
anticipated to play a substantial role. Such fluctuations 
are beyond the TDHF scheme. Not only TDHF cal- 
culations are, by construction, cylindrically symmetric 
around the beam axis. Such a limitation is not neces- 
sarily true in nature and the mean field dynamics could 
be not correct especially in presence of large fluctuations. 
Molecular dynamics contains all possible correlations and 
fluctuations due to the initial conditions (events). For the 
purpose of treating quantum- mechanical systems like tar- 



get atoms and molecules, we use classical equations of 
motion with constraints to satisfy the Heisenberg uncer- 
tainty principle and the Pauli exclusion principle for each 
event [g . In extending the study to the lower incident en- 
ergies, we would like to stress the connection between the 
motion of bound electrons and chaos. In fact, depending 
on the dynamics, the behavior of the electron(s) is un- 
stable and influence the relative motion of the projectile 
and the target. We could compare the D+d case to the 
gravitational 3-body problem, which has the same form of 
the equation of motion and it is nonintegrable |9j . For in- 
stance, the motion of asteroids around the sun perturbed 
by Jupiter becomes unstable, i.e., chaotic, depending on 
the ratio of the unperturbed frequencies of the asteroid 
and Jupiter. We discuss the enhancement factor of the 
laboratory cross section in connection with the integra- 
bility of the system by looking the inter- nuclear and elec- 
tronic oscillational motion. More specifically we analyze 
the frequency shift of the target electron due to the pro- 
jectile and the small oscillational motion induced by the 
electron to the relative motion between the target and the 
projectile. We show that the increase of chaoticity in the 
electron motion decreases the fusion probability. In this 
letter we will discuss the D+d case only because the sys- 
tem is fundamental to see its connection with chaos and 
has been well studied theoretically. We mention that the 
understanding of the fusion dynamics and fluctuations 
has a great potential for the enhancement of the fusion 
probability in plasmas for energy production. 

We denote the reaction cross section at incident energy 
in the center of mass E by <j(E) and the cross section 
obtained in absence of electrons by <jq{E). The enhance- 
ment factor f e is defined as 
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If the effect of the electrons is well represented by the 
constant shift U e of the potential barrier, following pilicj. 

(U e « E): 



fe ~ exp 
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where 17(E) is the Sommerfeld parameter [TT| . 
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We estimate the enhancement factor f e numerically 
using molecular dynamics approach; 
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where (r^p,) are the position, momentum of the par- 
ticle i at time t. £i = \J pfc 2 + mfc 4 , U(ri) and m, 
are its energy, Coulomb potential and mass, respectively. 
We set the starting point of the reaction at 10A inter- 
nuclear separation. Initially the electron is located in a 
Bohr orbit. To take into account the quantum mechani- 
cal feature of atoms, we put the constraints, i.e., Heisen- 
berg uncertainty principle and Pauli principle for atoms 
which have more than 2 bound electrons. It is performed 
numerically by checking Ar- Ap ~ h, for Heisenberg prin- 
ciple, and Ar • Ap ~ 27rft(3/47r) 2 / 3 , for Pauli blocking. 
Here Ar = |rj — r,-| and Ap = |pj — pj|. i and j refer to 
electrons and nuclei. More specifically, to get the atomic 
ground states, at every time step of the calculation, we 
calculate Ar • Ap for every pair of particles. If Ar ■ Ap is 
smaller (larger) than h, in the case of Heisenberg princi- 
ple, we change Yj and Pj slightly, so that Ar- Ap becomes 
larger(smaller) at the subsequent time step. We repeat 
this procedure for many time steps until no changes are 
seen in the energies and mean square radii of the atoms, 
similarly for the Pauli principle. The approach has been 
successfully applied to treat fermionic properties of the 
nucleons in nuclei and the quark system |8|. It can be 
extended easily in the case of the Heisenberg principle, 
as stated above. In this way we obtain many initial con- 
ditions which occupy different points in the phase space 
microscopically. Notice that in the D+d case that we 
investigate here, since one electron is involved, only the 
Heisenberg principle is enforced for each event. We ob- 
tain — 13.56 eV as the binding energy of deuterium atom 
and 0.5327 A as its mean square radius. These values can 
be compared with the experimental value of —13.59811 
eV [13 and Bohr radius R B = 0.529 A, respectively fl3|. 

In order to treat the tunneling process, we define the 
collective coordinates R co " and the collective momentum 
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where rp,rp (pt,Pp) are the coordinates (momenta) of 
the target and the projectile nuclei, respectively. When 
the collective momentum becomes zero, we switch on the 
collective force, which is determined by Fp " = p co11 and 
Fy W = —~p cM i to enter into imaginary time 0]. We 
follow the time evolution in the tunneling region using 
the equations, 
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where r is used for imaginary time to be distinguished 
from real time t. r^(p) and P^(p) are position and mo- 
mentum of the target (the projectile) during the tunnel- 



ing process respectively. Adding the collective force cor- 
responds to inverting the potential barrier which becomes 
attractive in the imaginary times. The penetrability of 
the barrier is given by [LH 



11(E) = {l + exp{2A{E)/h)Y 
where the action integral A(E) is 
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r a and r\, are the classical turning points. The internal 
classical turning point is determined using the sum of 
the radii of the target and projectile nuclei. Similarly 
from the simulation without electron, we obtain the pen- 
etrability of the bare Coulomb barrier Ho(E). 

Since nuclear reaction occurs with small impact param- 
eters on the atomic scale, we consider only head on col- 
lisions. The enhancement factor is thus given by eq.Q, 



f e =U(E)/Il (E) 



(8) 



for each event in our simulation. Thus we have an en- 
semble of f e values at each incident energy. 

In Fig. ^ the upper panel shows the incident energy 
dependence of the enhancement factor for the D+d reac- 
tion. The averaged enhancement factors f e over events 
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FIG. 1: Enhancement factor as a function of incident center- 
of-mass energy for the D+d reaction (upper panel). The 
corresponding E 2 // e (stars) and a power low fit(dashed line) 
(lower panel). 
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in our simulation are shown with stars and its variance 
£ = (/| — (/ e ) 2 ) with error bars. In the figure we 
show also several estimations of the enhancement fac- 
tor by the latest analysis of the experimental data using 
quadratic(dotted) and cubic (dot-dashed) polynomial fit- 
ting with the screening potentials U e — 8.7 eV and 7.3 
eV respectively. The dashed curve shows the enhance- 
ment factor in the adiabatic limit /i^" ' for an atomic 
deuterium target and it is obtained by assuming equally 
weighted linear combination of the lowest-energy gerade 
and ungerade wave function for the electron, reflecting 
the symmetry in the D+d, i.e., 

fe AD) = \ (exp (71-77 + exp(7r?7(i;)^)) , where 

U [ e 9) = 40.7 eV and ui u) = 0.0 eV H 0. In the low 
energy region the enhancement factor is more than 50. 
However the averaged enhancement factor does not ex- 
ceed the adiabatic limit. We performed also a fit of our 
data using eq. an d obtained U e = 15.9 ± 2.0 eV. This 
value, between the sudden and the adiabatic limit, is in 
good agreement with TDHF calculations [S0- 

The ratio E 2 / f e versus incident energy is plotted in the 
lower panel of Fig. ^ The numerical results(stars) dis- 
play a self similar behavior which is well fitted by a power 
law with exponent —2.6. In the high energy limit the 
ratio approaches zero, i.e., the f e distribution becomes a 
(5-function (S = 0) and f e — > 1: no effects due to the elec- 
tronic motion. In the low energy limit S 2 // e 3> 1, which 
implies a very sensitive dependence of the dynamics on 
the initial conditions, i.e., occurrence of chaos. Thus it 
is the motion of the electron which sensitively couples to 
the relative motion of the ions. 

Similar to the gravitational 3-body problem, we look 
at the oscillational motions of the particle's coordinates 
as the projection on the z-axis (the reaction axis). We 
denote the z-component of rr,rp and r e as zt,zp and 
z e , respectively. Practically, we examine the oscillational 
motion of the electron around the target zx e — z e — zt 
and the oscillational motion of the inter-nuclear motion, 
i.e., the motion between the target and the projectile, 
z s = zt + zp, which essentially would be zero due to the 
symmetry of the system, if there were no perturbations. 
In Fig. |21 these two values are shown for 2 events, which 
have the enhancement factor f e — 170.8 (ev. A), and 
f e = 6.5 (ev. B), at the incident energy E cm = 0.15 keV. 
The panels show the z Si zt c as a function of time. The 
stars indicate the time at which the system reaches the 
classical turning point. It is clear that in the case of ev. B 
the orbit of the electron is much distorted from the unper- 
turbed one than in ev. A. Characteristics of z s are that 
(1) its value often becomes zero, as it is expected in the 
un-perturbed system, and (2) the component of the devi- 
ation from zero shows periodical behavior. It is remark- 
able that the amplitude of the deviation becomes quite 
large at some points in the case of ev. B which shows the 
small enhancement factor. Note that in event B one ob- 
serves clear beats, i.e., resonances. Thus for two events, 
with the same macroscopic initial conditions, we have a 



completely different outcome, which is a definite proof of 
chaos in our 3-body system. We can understand these 
results in first approximation by considering the motion 
of the ions to be much slower than the rapidly oscillat- 
ing motion of the electrons. Thus we can consider the 
electron acting as an external force F e = Fq cos(u>nt + 6) , 
where Fq is the amplitude of the force and wh is the (hy- 
drogen) frequency. This will induce [15j a perturbation 
on z s ~ F e / (iuift, where [i is the ions reduced masses. No- 
tice how the amplitude of z s is actually reduced from the 
amplitude of zx e of a factor 1 ~ 10~ 4 , i.e., the ratio of the 
electron to the ion mass. We stress that this simple es- 
timation is more relevant for case A. In fact the motion 
of the electron is not decoupled from the inter-nuclear 
motion and from energy conservation we can expect that 
when \z s \ is maximum, \zx e \ is minimum as observed in 
Fig0 cas c B. Thus more generally, one should consider 
a perturbation where Fq — » Fo(t) and o->h —* wh(0- The 
time dependence of the perturbation leads, as it is well 
known, to the occurrence of chaos as for parametric res- 
onance From the Fig|21we can deduce the following 
important fact. If the motion of the electron is initially 
in the plane perpendicular to the reaction axis, the en- 
hancement factor is large, case A (notice \ztb\ Rb, i-e., 
the Bohr radius, at t ~ 0). On the other hand if there 
is a substantial projection of the electron motion, as in 
case B(thc amplitude of \zTe\ ~ Rb at t ~ 0), on the 
reaction axis the enhancement factor is relatively small 
because of the increase of chaoticity. This suggests that 
if one performs experiments at very low bombarding en- 
ergies with polarized targets, the enhancement factor can 
be controlled by changing the polarization. The largest 
enhancement with targets polarized perpendicularly to 
the beam direction. We notice in passing that event A 
is a case where cylindrically symmetry is approximately 
satisfied, because the electron motion is practically on 
the xy— plane. This case gives a U e =19.5eV closest to 
the adiabatic limit and to the TDHF result 00. 

In Fig.[3]we show the maximum amplitude of z s for 200 
events as a function of the enhancement factor f e at the 
incident energy E cm = 0.15 keV. Here we observe an evi- 
dent correlation between the two values: the events which 
give relatively large enhancement factors correspond to 
the small maximum amplitude of z s , and to the contrary, 
the events with relatively small enhancement factors indi- 
cate large maximum amplitudes. Since the motion of the 
ions is coupled to that of the electrons, it implies that an 
amount of energy is transferred from the relative motion 
of the ions to the electron thus reducing the probability 
of fusion. On the contrary small amplitudes imply that 
the electronic configuration is closer to the one in the g.s. 
of the compound system. Thus the binding energy of the 
electron is converted into the relative energy of the ions 
increasing the fusion probabilities. 

In conclusion, we discussed the penetrability of the 
Coulomb barrier by using molecular dynamics simula- 
tions with constraints and imaginary time. We have 
shown that both the enhancement factor and its vari- 
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FIG. 2: The oscillational motion of the electron around the target (lower panels) and the inter-nuclear motion (upper panels) 
as a function of time, in atomic unit, for two events, with large / e (ev. A) and small / e (ev. B), for the D+d reaction at the 
incident energy 0.15keV. The inter-nuclear separation is 10A at t = 0. 




ance increase as the incident energy becomes lower. How- 
ever we obtained the averaged screening potential smaller 
than the value in the adiabatic limit, while from fluctua- 
tions some events clearly exceed such a limit. We pointed 
out that there is an evident correlation between the os- 
cillational motion of the inter-nuclear separation and the 
magnitude of the enhancement factor of the cross sec- 
tion. The chaoticity of the electron motion affects the 
enhancement factor of the cross section. We suggest to 
perform experiments on fusion at very low energies with 
polarized targets. 



FIG. 3: The maximum amplitude of z s for 200 events as a We acknowledge valuable discussions and suggestions 

function of the enhancement factor f e for the D+d reaction witn p rofs D M Brink, G. Fiorentini, C. Rolfs, C. Spi- 
at the incident energy 0.15 keV. taleri and N _ Takigawa . 
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